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ABSTRACT 

Recent experiences in utilizing the global optimization 
methodology, based on polynomial and neural network 
techniques, for fluid machinery design are summarized. Global 
optimization methods can utilize the information collected from 
various sources and by different tools. These methods offer 
multi-criterion optimization, handle the existence of multiple 
design points and trade-offs via insight into the entire design 
space, can easily perform tasks in parallel, and arc often 
effective in filtering the noise intrinsic to numerical and 
experimental data. Another advantage is that these methods do 
not need to calculate the sensitivity of each design variable 
localls However, a successful application of the global 
optimization method needs to address issues related to data 
requirements with an increase in the number of design variables, 
and methods for predicting the model performance. Fxamples of 
applications, selected from rocket propulsion components, 
including a supersonic turbine and an injector element, and a 
turbulent flow diffuser are used to illustrate the usefulness of the 
aloha! optimization method. 

t. INTRODUCTION <& OV ERVIEW 

Modern computational and experimental fluid 
dynamics tools have matured to a stage where they can provide 
substantial insight into engineering processes involving fluid 
Hows. 1 his can help analyze the fluid physics as well as improve 
i he design ol practical devices. In particular, rapid and 
continuous development in the technology of fluid machinery 
demands that new design concepts be regularly proposed to 
meet goals For increased performance, robustness and safety 
while concurrently decreasing cost. These stringent goals are 
forcing consideration of design variables over ranges and in 
combinations not typically employed, thereby increasing the 
design complexity. Objective and efficient evaluation of these 
new and complex designs can be facilitated by development and 
implementation of systematic optimization methods. To date. 
The majority of the effort in design optimization of fluid 
dynamics has relied on gradient-based search algorithms (Baysal 
and FJeshaky 1992, Lambert et al. 1995, Reuther et ai. 1999). 
These methods work iteratively through a sequence of local sub- 
problems. which approximate objective and constraint functions 
for a sub-region of the design space. e,g., by linearization using 


computed sensitivities. Major challenges for these optimization 
approaches are the robust and speedy computation of sensitivity 
coefficients (Elbanna and Carlson 1994. Dadonc et al. 2000). 

Polynomial-Based Response Surface 

Recently we have witnessed increasing interests in 
utilizing the global optimization methods, which avoids the need 
to compute design sensitivities altogether. In this regard, the 
response surface methodology (RSM) has gained attention 
because minimum eflbrt is required lor software interfacing, 
which facilitates the integration of analysis and design codes. 
RSM is particularly suitable for subsystem approximation within 
multidisciplinary' design optimization (Sobieszczanski-Sobieski 
and Malika 1997) because it can utilize information collected 
from various sources and by different tools. It can also offer 
multi-criterion optimization, handle the existence of multiple 
design selections and related trade-offs, and address the noises 
intrinsic to numerical and experimental data. A main advantage of 
RSM lies in its robustness and intelligibility. Robustness and the 
smoothness of the approximations in noisy environments are 
achieved by performing extra analyses compared to the number 
of regression coefficients. This is a distinct advantage over 
derivative-based search algorithms, which may encounter 
difficulties in the presence of spurious local optima. This 
generality allows the consideration of information at varying 
levels of breadth ( Le.. number of design \ariables) and depth (i.c.. 
details of the design space). The significance of the individual 
design parameters can also be assessed directly by the global 
model. Another feature favoring RSM is its suitability for parallel 
computation. 

The polynomial-based response surfaces are commonly 
employed in global optimization. An n' h order polynomial in each 
dependent variable is generated by standard least-squares 
regression (Myers and Montgomery 1995). Polynomial 
coefficients are obtained from a linear regression scheme. Tn 
addition to polynomials, neural networks iNN) can also be used 
to construct a response surface model. With a polynomial 
approach, there is a cost associated with RSM. The computational 
demands grow fast with increasing number of design variables, 
which is termed as the ** curse of dimensionality" . For example, to 
construct a second-order polynomial of A’ design variables, the 
number of coefficients to be fixed are ( S - l)(N+2)/(2!). A cubic 
model will require (N+l)(N+2)(N+3) (3!) coefficients. In 
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addition, the predictive capability of RSM is greatly influenced 
by the distribution of the sampling points in the design space 
(Unai el al. 1997, 1998). 

Statistical techniques are also available for identifying 
polynomial coefficients that are not well characterized by the 
data. For example, a stepwise regression procedure based on t- 
statistics is often used to discard terms and improve the 
prediction accuracy. The t-statistic, or l-ratio. of a particular 
coefficient is given by the value of the coefficient divided by the 
standard error of the coefficient. The quality of fit of different 
surfaces can be evaluated by comparing the adjusted root mean 
square error. a it , defined as: 



Mere e t is the error at i ,h point of the training data, n is the 
number of training data points and n (i is the number of 
coefficients, file measure of error given by a, is normalized to 
iccount for the degrees of freedom in the model. Thus a. 
accounts for the nominal effect of higher order terms providing a 
better overall comparison among the different surface fits. 

The accuracy of the models in representing the design 
<pacc is further gauged by comparing the values of the objective 
function at test design points, different from those used to 
generate the fit. with the empirical solution. The root mean 
'Muure error, a : for the test set is gi\en by: 



in this equation is I he error at the lest point and m is the 
anal number of test points. 

Fvpicaily. the entire design data are divided into two 
groups. The major portion, called the training set. is used to llx 
the global model, and the minor portion, called the test set. is 
used to lest the fidelity of the model. In short, the variation 
between the response surface (RS) and the training data, as 
given by Fq.(l), is normally used to judge the performance of 
die tit. As indicated in Hq. (2), additional test data can be 
employed to evaluate the performance of different polynomials 
■ er design points not used during the training phase 

in fluid dvnamics applications, the RSM has been 
applied to a growing number of cases, including high speed civil 
transport iKnill cl ai. 1999), airfoil shape optimization (Rai and 
Madavnn 1998. 2000. Madavnn et al. 1999). diffuser shape 
optimization (Madsen et al. 2000). preliminary design of 
supersonic turbine, ( Pap i la et al. 2000). and injectors (Shyy et al 
1999. Tucker etai. 1998). 

Data Selection in Design Space 

In order to reduce the size of the data needed for 
constructing the global model, it is important to apply a sound 
strategy for selecting design points (Malika el al. 1998). This 
strategy is referred to as the design of experiments. For example, 
face centered composite design (FCCD) is a popular approach 
or selecting the training data to aid the construction of a 
response suriace model, fin tv three-dimensional space. FCCD 
creates a design space composed of eight corners of the cube. 


one at the center of each of the six faces and one at the center of 
the cube. Therefore, this yields (2 ? *+2N+I) points, where N is the 
number of design variables. Due to this scaling rule, FCCD is 
more effective when the number of design variables is modest, 
say, no larger than 5 or 6, but not a good choice for problems with 
large number of design variables. 

The orthogonal array (OA) is a fractional factorial 
matrix that assures a balanced comparison of levels of any factor 
or interaction of factors. Because the points are not necessarily at 
vertices, the analytical tools can be more robust using the 
orthogonal array. The OA can significantly reduce the number of 
data required to construct a global model. To further reduce the 
data size, one can choose to rank the different point selections 
according to D-optimality criterion. This approach minimizes the 
generalized variance of the estimates, and can reduce the 
sensitivity of the response surface with respect to noise. In Papila 
et al. (2000), alternative representations of the design space are 
performed for supersonic rocket turbine design by using FCCD 
coupled with the D-Optimal treatment and OA. FCCD produces 
"-data for the single-stage turbine with six design variables. 
With 1 1 design variables for the two-stage turbine design. FCCD 
fields 2.071-data. For three-stage turbine having 15 design 
variables. FCCD generates 32.799-data based on the formula of 
2 a+ 2.Y+/, demonstrating the "curse of dimensionality M . For such 
cases methods like OA can be applied to reduce the number of 
data in an efficient way. 

To enhance the effectiveness of the RSM. as 
demonstrated by Papila cl al. (2000). a multi-level approach is 
very useful. One can construct a global response surface and 
■ Jcntily regions expected to be of favorable performance. Then, 
ut these selected regions, one can retine the global model and 
identify the optimal design points with substantially higher 
tecuraev . 

\euraf \etworhs 

As already mentioned. NN techniques have also been 
used to generate surrogate models representing data obtained 
from simulations based on complex numerical and experimental 
schemes. Relevant papers by Carpenter & Barthclemy (1993). 
Nikolaidis el al. ( 1998). Grcenman and Roth ( 1998). Madavan et 
al. (1998). Rai and Madavan (1998. 2000). Shyy et al. (1999), 
Papila et al. (1999) and Vaidyanalhan et al. (2000) include 
comparative studies of polynomials and neural networks for data 
handling. NN is highly flexible in functional form and hence can 
oi'er 'diimlicant potential for complex functions that cannot be 
adequate!} approximated by polynomials. NN can be effectively 
used in two ways. First, it can be used in conjunction with RSM. 
For example, in complex regions of the surface. NN can he 
trained using the existing data and then can be used to generate 
additional data thereby enhancing the available information in 
that particular area. Shyy et al. (1999) have demonstrated that this 
approach can help improve the fidelity of the polynomial-based 
response surface model. Alternatively, NN can generate data that 
can be directly used by the optimizer. With either a polynomial- 
based response surface method or a neural network method, the 
design optimization is conducted h\ first constructing a global 
representation of the design space, followed by a search for 
optimal designs. 

While back-propagation neural networks (BPNN) have 
been largely employed in reported studies for NN -based design 
optimization, the radial-basis neural networks (RBNN) have 
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advantages lhat deserve to be investigated more closely. RBNN 
are multi-layer networks with hidden layers of radial basis 
transfer function and a linear output layer. The training of the 
network is a cyclic process and the weights and biases of the 
nodes of the network are adjusted until an accurate mapping is 
obtained. RBNN may require more neurons than BPNN, but 
they can be designed faster than the latter. Various parameters 
need to be evaluated to design an RBNN. Tor example, a spread 
constant is needed as a design parameter; i.e., the radius of the 
basis in the input space to which each neuron responds. One can 
design a network with zero error on the training vectors by 
generating as many radial basis neurons in the network as there 
arc input vectors. A more compact design in terms of network 
size is obtained by generating one neuron at a time to minimize 
the number of neurons required. At each cycle/epoch, a neuron 
is added to the network until a user specified error goal is 
achieved or until the network has generated the maximum 
number of neurons possible. In such a case, an error goal and 
the spread constant need to be specified. The lest data helps to 
evaluate the accuracy of the networks with varying spread 
constants. Thus, the test data are part of the evaluation process 
and help in selecting the best NN. The NN technique can 
encounter problems due to the nature of the data, as well as the 
construction of the neuron characteristics. If the training dataset 
is noisy, and the network is trained without proper filtering 
features, false optima due to over-lilting will occur. Raise 
optima can also be introduced when inadequate training data set 
*s used. 

In this paper, we summarize our recent experiences in 
utilizing the global optimization methodology related to 
pol\ nomial-based RSM and NN. Then, we highlight three 
physical examples: a supersonic turbine and an injector element 
for rocket propulsion applications, and shape optimization of a 
turbulent How diffuser. To reach a successful optimal design, 
one often needs to consider the issues related to (i) selection of 
appropriate representation of the design space for constructing 
the global model, ( ii) employment of the statistical and testing 
tools to assess appropriate global models, (iii) multi-criterion 
optimization. < iv) existence of multiple design selections and 
related trade-offs, and (v) consideration of noises intrinsic to 
numerical and experimental data. There is no space to address 
all those matters in this paper, but detailed information can be 
found in the references cited. 

:. PHYSICAL APPLICATIONS 


the chamber wall heat flux, Q. Shyy et al. ( 1999) and Tucker ct al 
(1999) used a geometric mean to combine these two objectives 
and maximized the resulting composite objective function, D. 

D=((d m y{d Q ]f o) 

w'here the normalized function such as ERE , takes the form: 


d /.///. — 


( ERE- A 
l B-A 


(4) 


where B is the target value and A is the lowest acceptable value. 
Figure 1 illustrates the roles of 5 and t in the desirability function 
for the case of maximizing a response. The desirabilities with 
s<<! imply that a product need not be close to the response 
target value, £, to he quite acceptable. However, a large value ofs 
implies that the product is nearly unacceptable unless the 
response is dose to B. 

Hie other way of expressing the composite function is 
to use a weighted sum of the objective function. The composite 
desirability function can then be expressed as 


l) = aj\ +aj 2 + •■■ (5) 

where D is the composite objective function and fs arc the non- 
normalized objective functions, flic arc dimensional 
parameters that control the importance of each objective function. 
In the case i>T Llie supersonic turbine a weighted sum of the two 
objectives, the efficiency, and the weight. W. has been used. 
The expression, in the context of the turbine gives the incremental 
’value of the payload. Apaw with the change in /Land //. The goal 
is to maximize J pay. which in turn results in maximizing ;; and 
minimizing /Land lienee the payload is maximized. 

D - ,:J pay - Ci( /]- J]/.) x\ ()()-( }\ - 1 F,,y ( 6) 

where }]= the calculated efficiency 
ijh = the baseline efficiency 
JF = calculated weight 
IF* = the baseline weight 

c t = the amount of payload increment capacity for any 
eifieienev earn 


Object i ve Functions 

When attempting to optimize two or more different 
objectiv e functions, conflicts between them arise because of the 
different relationships they have with the design variables. To 
solve this problem, a multi-objective approach is often 
employed. Here, competing objective functions are condensed to 
a single function by generating a composite objective function. 
The maximization of the composite function effectively 
optimizes the individual functions. The use of a response surface 
type of global model makes it straightforward to do such a 
multi-criterion optimization. W ithout response surfaces, it would 
have been a highly challenging task. 

To handle such a multi-criterion optimization task, an 
average of some form is normally used to represent the 
composite function, l or example, for the injector, the goal is to 
maximize the energy release efficiency, ERE while minimizing 


The baseline efficiency and weight are obtained with 
the existing design knowledge without benefiting from an 
optimization strategy. The goal of the optimization is to identify 
turbine configurations capable of delivering higher payload. 

Turbulent Flow Diffuser 

Figure 2 highlights the use of a response surface 
approximation for the optimum shape of a two-dimensional 
diffuser. The goal was to accomplish maximum pressure recovery 
by optimizing the wall contours. The How is incompressible and 
fully turbulent with a Reynolds number of 10', based on the inlet 
throat half-width, D. The overall geometry is defined by the ratio 
of inlet and outlet areas, and the diffuser length to height ratio, In 
-t-his study the length to height ratio is fixed at 3.0, and the area 
ratio at 2.0. The shape of the diffuser wall is designed for 


3 


optimum performance, with live design variables represented by 
B-splines. The CFD model is based on the full Reynolds- 
averaged Navier-Stokes equations, with the k-e two-equation 
turbulence model in closure form. At the inlet of the flow 
domain, a uniform flow distribution is specified. Detailed 
discussion of this study can be found in Madsen et al. (2000). As 
illustrated in Fig. 1, within the fidelity of the analysis tool, there 
are often multiple design points that meet the design objectives. 
It is interesting to note that different diffuser shapes are found to 
yield essentially the same performance. The response surface 
model is ideally suitable for such situations. 

Supersonic Turbine 

Next, we summarize our recent efforts in optimizing 
the preliminary design of a supersonic turbine suitable for a 
reusable launching vehicle (RLV) propulsion svstem Single-, 
two- and three-stage turbines arc considered with the number of 
design variables increasing in accordance with the number of 
stages, ( here arc 2 types of design variables: 

ii come trie input 

• mean diameter (/>) 

• last rotor annulus area (A tum ) 

• blade height ratio between the T 1 Mine and the last rotor 

blade 

• vane and blade axial chords 
Performance f 'anahics 

• ar\f 

• number of stages 

• blade row reaction 

• work split (if more than I stage is investigated) 

In the work presented by fiapila et al. (2000). there are 
o. i I. and 15 design variables lor the single-, two- and three- 
-tage turbines, respectively. In addition, there are 2 structural 
constraints, the blade centrifugal stress and the disk stress. The 
blade centrifugal stress was constrained by a limit placed on the 
product of the blade exit annulus area and the RP\i~ f. l.V / . I'he 
disk stress was constrained by a limit placed on the pitchline 
velocity (the product of the RPM and the mean radius). All the 
design variables involved in the design process are normalized 
by their respective baseline values. 

lor rocket engine applications, maximizing the 
chicle payload for a given turbine uperaiing condition is the 
ultimate objective. Any gain in turbine efficiency will he 
re lice led m a reduced propellant consumption, thus an increase 
m payload. However, higher turbine performance usually emails 
multistage designs, which arc heavier. To ascertain required 
predictive capability of the RSM, a two-level domain refinement 
approach has been adopted in the course of optimization (Papiia 
et al. 2000). First, a response surface was constructed for the 
enure design space. Since the accuracy of the response surface 
was less than satisfactory, a domain refinement was then 
adopted based on the initial optimization. 

An inspection of the optimal designs indicates that the 
two-stage turbine gives the best payload performance for this 
application. As the number of the stage increases, we see that 
efficiency improves whilq the weight increases also. According 
to the formula for Apay , the improvement in efficiency from 


two- to three-stage cannot compensate the penalty from higher 
weight. As shown in Figure 3, the mean diameter, speed, and the 
exit blade area exhibit distinct trends. Specifically, the diameter 
decreases, speed increases, and annulus area decreases with 
increasing number of stages, ft is interesting to observe that none 
of the selected values of the design parameters are toward the 
limiting values assigned, indicating that the optimal designs result 
from compromises between competing parametric trends. For 
such cases, a formal optimizer such as the present response 
surface method is very useful. The results indicate that the 
efficiency rises quickly from single stage to two stages but the 
increase is much less pronounced with three stages. A single- 
stage turbine performs poorly under the engine balance boundary- 
condition possibly due to a significant portion of fluid kinetic 
energy being lost at the turbine discharge of the single-stage 
design due to high stage pressure ratio and the high-energy of the 
working fluid, which is predominantly hydrogen. Figure 4 
summarizes filting/training and testing results of RBNN and 
polynomial-based approximations for Apay of a 2-stage turbine. 
1 he training data arc .selected hased on the orthogonal array 
approach. There are l I design variables. 249 training data (OA). 
and 78 testing data in both original and refined design spaces, 
(/entered on the optimal design point predicted using the 
information gathered for the entire design space, the reduced 
design space has 20% of the range original ly assigned to each 
design variable. In all plots, a perfect fit will result in a 45-degree 
line. 

I he effectiveness of the muiti-Ievel RSM approach can 
be observed in this figure by comparing the original and refined 
design space plots. W hile RBNN is obviously more accurate for 
the training data due to more number of adjustable parameters 
associated with the number of neurons, its predicting capability, 
as retlected in the testing data plots. mn\ not be as superior. 
Nevertheless. RBNN is attractive because it has the llcxibilitv to 
handle complex characteristics of a response surface and is 
straightforward to train because the computation is based on the 
linear regression analysis. 

Swirl Coaxial injector Element 

The injector considered is tor die combustor of an 
advanced liquid rocket engine. Different injector concepts have 
been proposed for this purpose, including shear coaxial type, 
impinging type, and swirl type {Shyv et al. 1999, Tucker et al. 
1999). Hie swirl coaxial element lias been used somewhat 
■ paringly ;n the U.S.A.. but has been wide!) used in Russia 
because of its reported ability to perform well over a large throttle 
; ange ((till and Nurick 1976). Figure 5 shows that the GO : , 
Rowing in the center post of the element, exits from the element 
with both radial and axial velocity components. This effect is 
achieved by introducing the G0 2 tangentially into the center post 
through small slots. The empirical design methodology of 
Calhoon et al. (1973) uses the oxidizer pressure drop, AP () , fuel 
pressure drop. AP h combustor length. / Vft/nA , and the full cone 
swirl angle. 6T as independent variables. Due to stability 
considerations for this injector design, the AP W range is set to 10- 
20% of the chamber pressure, while the Afj range is set to 2-20% 
of chamber pressure. The combustor length, defined as the 
distance from the injector to the end of the barrel portion of the 
chamber ranges from 2-8 inches. The lull cone swirl angle is 
'allowed to vary from 30-90°. The dependent variables modeled 
are ERE (energy release efficiency, a measure of element 
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performance), wall heat flux, injector heat flux. O m! , relative 
combustor weight, W n , h and relative injector cost, C rt . f . 

With the multi-criterion optimization, one can assess 
the effect of certain aspects of the design during the optimization 
process. The set of results shown in Table 1 facilitate the 
illustration. The baseline results in Table 1 are obtained by 
placing equal weight for all design variables. Alternatively, in 
Case I, we place more emphasis on reducing the wall and 
injector face heat fluxes. Desirability functions for both of these 
variables are given increased weights (5 and 10, respectively). 
Since lower heat fluxes tend to increase component life, 
weighting these two variables is equivalent to emphasizing a 
life-type issue in the design. Since Q„ is already at its minimum 
value, it remains fixed. As expected. & is decreased which 
decreases the value of 0 HI/ by almost 35%. The lower value of <9 
also produces a lower ERE. Both propellant pressure drops and 
the combustor length are increased to mitigate the drop in ERE 
The increases in L atm h and AP, cause an increase in Jf' w and 
respectively. The emphasis on life extracts the expected penally 
on performance. Additionally, for the current model, there are 
also flight weight and cost penalties. 

fhe results for Case 2 are obtained b> emphasizing 
maximization o \ ERE and minimization of lf r /vi vv ilh desirabilitv 
weightings of 10 and 5. respectively. Increased weighting for 
ihese two variables is equivalent to emphasizing a thrust to 
weight goal for the injcctor/ehamber. The relative chamber 
length is shortened to slightly lower ERE is maximized by 
increasing the Ci(T swirl angle by a factor of almost 2.5 and 
increasing AP, b\ over 35 %. The value of ERE increases by 
over one percent As noted earlier, increasing rt leads to 
increased injector heal tlux. For this case, emphasis on thrust 
and weight tends to have an adverse affect on (J,„ r Relative cost. 

( . lor the current model, is also increased significantly. Figure 

o offers a visual illustration of the trends exhibited in various 
scenarios between a swirl type and an impinging t\pe injector. 
While the qualitative features between the two injectors are 
consistent, the present optimization model allows the inspection 
of the parametric variations direct!}. 

3. CONCLUSIONS 

The global optimization method can be particularly 
effective with either a polynomial-based response surface or a 
neural network when information from different computational, 
experimental. and analytical sources needs to be assembled. In 
•m> paper, we summarize our recent experiences in utilizing the 
optimization methodology for optimizing tasks related to the 
preliminary design of a supersonic turbine, the design of an 
injector element, and the shape optimization of a turbulent flow 
diffuser. A successful optimal design often needs to address the 
issues related to the selection of appropriate training data for 
constructing the global model, employment of the statistical and 
testing tools to identify appropriate global models, multi- 
criterion optimization, existence of multiple design selections 
and related trade-offs, and consideration of noises intrinsic to 
numerical and experimental data. 

4. ACKNOWLEDGMENTS 

File work presented has been supported b\ NASA 
Marshall Space Flight CenteiN The paper is a summary of our 


collaboration with our colleagues, as cited in the references given 
below. 

5. REFERENCES 

Baysal, O. and Eleshaky, M.E., 1992 "Aerodynamic 
Design Optimization Using Sensitivity Analysis and CFD", AlAA 
Journal Vol.30 pp.7 1 8-725. 

Calhoon, D., Ito, J. and Kors, D., 1973 “Investigation of 
Gaseous Propellant Combustion and Associated Injector- 
Chamber Design Guidelines," Aerojet Liquid Rocket Company, 
XASA CR-l 21234, Contract NAS3-i33 7 9 . 

Carpenter, W.C. and Barthelemy, J.-F.M, 1993 “A 
comparison of Polynomial Approximations and Artificial Neural 
Nets as Response Surface," Structural Optimization , Vol. 5, pp. 
166-174. 

Dadone, A., Valorani, M. and Grossman. B., 2000 
"Smoothed Sensitivity Equation Method for Fluid Dynamic 
Design Problems," AIAA Journal , Vol. 38, pp. 418-426. 

Elbanna. II.M. and Carlson, L.A.. 1994 "Aerodynamic 
Sensitivity Coefficients Using the Three-Dimensional Full 
Potential Fquations", J. Aircraft Vol.3l. pp. 1071-1077. 

Gill. G. S.. and Nurick. W. 11.. 1976 Liquid Rocket 
Engine Injectors , NASA SP-8089. 

Greenman, R. M. and Roth. K.R.. 1998 "High-Lift 
Optimization Design Using Neural Network." Proceedings of 
DETC98, AS ME Computers in Engineering Conference . Atlanta. 
Georgia, DETC98/CIE-6006. 

Greenman. R. M, and Roth. K.R.. 1999 "Minimizing 
Computational Data Requirements For Multi-Element Airfoils 
1 sing Neural Networks." 3 m 'th AIAA Aerospace Sciences Meeting 
and Exhibit . Reno, Nevada, AIAA 99-0258. 

Malika. R.T.. Scott. F,P and Cruz. J R.. 1998 

"Optimization and Fxpcrimenis: A Survey." Applied Mechanics 
Review, Vol. 51. pp. 435-448. 

Kn i 11. D.L.. Giunta, A. A., Baker. C.A., Grossman, B.. 
Mason. W.II., I laftka. R.T. and Watson. FT.. 1999 "Response 
Surface Models Combining Linear and Euler Aerodynamics for 
Supersonic Transport Design," Journal of Aircraft, Vol. 36. pp. 
75-86. 

Lambert. IF A., Lccordix, J.L., and Braibant, V. 1995 
"Constrained Optimization of Nacelle Shapes in Euler Flow using 
Semi-analvtical Sensitivity .Analysis" Structural Optimization . 
Vol. 10. pp. 239-246. 

Madavnn. N.K. Rai. M M. and Huber. F W.. 1999 
"Neural Net-Based Redesign of Transonic Turbines lor Improved 
I ’nsteady Aerodynamic Performance." aHA/SAE/ASME/ASEE 
3 5 th Joint Propulsion Conference. AIAA 99-2522. 

Madsen, .1.1., Shyy, W. and I laftka, R.T. 2000 
"Response Surface Techniques for Diffuser Shape Optimization ," 
accepted for publication in AIAA Journal. 

Myers, R.H. and Montgomery 7 , D.C., 1995 Response 
Surface Methodology (Pro-cess and Product Optimization Using 
Designed Experiments). Wiley, New York. 

Nikolaidis. E.. Long, L., and Ling. Q., 1998 "Neural 
Networks and Response Surface Polynomials for Design of 
Vehicle Joints," 39th AIAA/ AS ME ASCE AHS/ASC Structures. 
Structural Dynamics, and Materials Conference , Long Beach. 
California, AIAA-98-I777. 

Papila. N.. Shyy. W.. Fiiz-Co>. N. and Haftka, R.T, 
-t9 9 9 "Assessment of Neural Net and Polynomial-Based 


Techniques for Aerodynamic Applications,” 17th Applied 
Aerodynamics Conference, Norfolk, VA, AIAA 99-3 167. 

Papila, N., Shyy, W., Griffin, L.„ Huber. P. and Tran, 
K., 2000 ’'Preliminary Design Optimization for a Supersonic 
Turbine for Rocket Propulsion,” A I A A/S A E/A SME/A SEE 35 th 
Joint Propulsion Conference , Paper No. 2000-3242. 

Rai, M.M. and Madavan, N.K, 1998 "Aerodynamic 
Design Using Neural Networks,” y th AIAA/USAF/ 
NASA/ISSMO Symposium on Multidisciplinary Analysis and 
Optimization , St. Louis, AIAA 98-4928. 

Rai, M.M. and Madavan, N.K, 2000 "Improving the 
Unsteady Aerodynamic Performance of Transonic Turbines 
Using Neural Networks," 38th AIAA Aerospace Sciences 
Meeting and Exhibit , Reno, Nevada, AIAA 2000-0169. 

Reuther, I, Alonso, J.J., Rimlingcr, M.J. and Jameson, 
A.. 1999 "Aerodynamic Shape Optimization of Supersonic 
Aircraft Configurations via an Adjoint Formulation on 
Distributed Memory Parallel Computers.” Computers and 
Fluids , Vol. 28. pp. 675-700. 

Shyy. W.. Tucker. ILK. and Vaidyanathan. R., 1999 
"Response Surface and Neural Network techniques for Rocket 
Fngine Injector Optimization.” AlAA.SAE/ASME.lSEE 35th 
Joint Propulsion Conference . AIAA 99-2455. 

Sobieszczanski-Sobieski. J. and Ilnftku. R.T., 1997 
"Multidisciplinary aerospace design optimization: survey of 
recent developments". Structural Optimization . Vol. 14. pp. I- 
23. 


Tucker, P.K., Shyy, W. and Sloan, J.G., 1998 “An 
Integrated Design/Optimization Methodology for Rocket Engine 
Injectors,” AlAA/SAE/ASME/ASEE 34th Joint Propulsion 
Conference , AIAA-98-3513. 

Tucker, P. K., Shyy, W. and Vaidyanathan, R., 1999 
"Optimization of a G0 2 /GH 2 Impinging Injector Element,” NASA 
l(f Thermal and Fluids Analysis Workshop , September 13-17, 
Huntsville, AL. 

Tucker, P. K., Shyy, W. and Vaidyanathan, R., 1999 
"Optimization of a G0 2 /GH 2 Swirl Coaxial Injector Element,” 
Penn State University Propulsion Engineering Research Center 
/ i !h Annual Symposium on Propulsion . November ( 1 999). 

Unal, R., Braun, R. D., Moore, A.A., and Lepsch, R.A., 
1997 "Response Surface Model Building Using Orthogonal 
Arrays for Computer Experiments," 19 th Annual International 
Conference of the International Society of Parametric Analysis , 
New Orleans, Louisiana. 

Unal, R., Lepsch. R. A., and McMiliin, M. L, 1998 
"Response Surface Model Building and Multidisciplinary' 
Optimization using D-Optimal Designs.” 

/ lAA/CSA F/NASA/ISSMO Symposium on Multidisciplinary' 
I n a lysis and Optimization, Paper No. 98-4759. 

Vaidyanathan. R.. Papila. N . Slny. W.. Tucker. P.K.. 
Griffin, L. W., Haftka. R.T, and Fitz-Coy. N.. "Neural Network- 
Based and Response Surface-Based Optimization Strategies for 
Rocket Fngine Component Design.” 8th 
1 1.4 A/USA F/NASA/ISSMO Symposium <)n \ lultidisciplinary 

Analysis and Optimization. Paper No. 2000-4880 


Table 1 Effect of emphasizing life and performance issues on optima! injector design. 
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Figure 2. C ontour plot of response surface for diffuser design. Solid circle indicates the optimal region, The hatched part of the feasible 
-pace comprises designs with performance within \% of the optimal. Corresponding shapes are indicated to the right. Flic results 

indicate that multiple design points meet the goal. 
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Figure 3. Effect of the number of turbine stage on optimum output parameters; ?], IV, and Apav and optimum design parameters; D, 
RPM. and A omi calculated for Jpav-based optimization. All design variables are normalized by the baseline values. 
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Figure 4. Comparison of NN and polynomial-based representations for 2-stage supersonic turbine. Plotted are the original and predicted 
\ alucs of Apciy% A perfect tit will result in a 45-degree line. The training and testing data are selected based on the orthogonal arrays with 
i he D-opiimat criterion v fhere are I I design variables. 249 training data (OA). and 78 testing data in both original and refined design 

spac.cs. (The values for dpay are normalized by the baseline value) 
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